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A single vortex overcoming the surface barrier in a mesoscopic superconductor with lateral dimensions of 
several coherence lengths and thickness of several nanometers provides an ideal platform to study thermal ac- 
tivation of a single vortex. In the presence of thermal fluctuations, there is non-zero probability for vortex 
penetration into or exclusion from the superconductor even when the surface barrier does not vanish. We con- 
sider the thermal activation of a single vortex in a mesoscopic superconducting disk of circular shape. To obtain 
statistics for the penetration and exclusion magnetic fields, slow and periodic magnetic fields are applied to the 
superconductor. We calculate the distribution of the penetration and exclusion fields from the thermal activation 
rate. This distribution can also be measured experimentally, which allows for a quantitative comparison. 

PACS numbers: 74.78.-w, 74.78.Na 



I. INTRODUCTION 

Thermal activation of macroscopic quantum objects over an 
energy barrier is ubiquitous in condensed matter physics. One 
well established example is the thermally assisted escape of 
phase particle in a Josephson junction when the bias current 
is close to the critical current corresponding to a vanishing 
barrier[ 1]. Combined with analytical calculations and exper- 
iments, a quantitative description has been achieved in this 
case. Magnetic vortices as topological excitations in super- 
conductors are difficult to excite thermally because of their 
large self-energy compared to that of phase particle in Joseph- 
son junctions. For a large of number of vortices the probabil- 
ity of thermal activation increases, which allows for an ex- 
perimental observation. Thermal activation of bundles of vor- 
tices, such as creep motion of vortices[2|, has been observed 
long time ago. However, no detailed quantitative experimen- 
tal study of thermal activation of a single vortex has been re- 
ported so far because the probability of thermal activation for 
a single vortex is usually small. 

Vortices overcoming the surface barrier in finite-size su- 
perconductors has attracted considerable interest, because of 
the possible observation of thermal activation of a quantum 
object. For a superconductor with finite size in lateral di- 
mensions, there exists a surface barrier known as the Bean- 
Livingston (BL) barrierQ. When the applied magnetic field 
increases, the BL barrier diminishes and finally disappears at 
the penetration field H p , where vortex enters the supercon- 
ductor. The penetration field is larger than H c \, the ther- 
modynamic lower critical field for type II superconductors 
without surface. For the same reason, the exclusion field 
H e , when vortex leaves the superconductor, is smaller than 
H c \. Thermally assisted penetration of vortices in a thick su- 
perconductor was investigated theoretically by Petukhov and 
Chechetkin decades ago. [4] They found that the thermal acti- 
vation is practically impossible for vortex. For high-r c super- 
conductors, Kopylov et. al. found experimentally that the 
penetration field decreases when the sample temperature is 
increased[5 1. The thermal effect on the penetration and exclu- 
sion field in high-r c superconductors has also been studied in 



Refs. [6-8 1. The importance of thermal fluctuations in these 
cases can be estimated using the Ginzburg number G 3D = 
\( k B T c /[H c (Q)^(Q) } ]) 2 , where k B is the Boltzmann constant, 
T c the critical temperature, H c {0) the thermodynamic critical 
field and £(0) the coherence length at zero temperature[9|. 
For conventional superconductors, G 3 ° ~ 1CT 10 and ther- 
mal fluctuations are weak; while for high-r c superconductors, 
G 3D » 1CT 2 and thermal fluctuations become important. 

For a large size superconductor studied in Refs. EHH), 
many vortices penetrate into the superconductor simultane- 
ously once the applied field reaches the penetration field. To 
observe single vortex penetration, the lateral size of the su- 
perconductor should be of several coherence lengths only. In 
these mesoscopic superconductors, a single vortex enters or 
exists when the applied magnetic field is tuned. Due to the 
advances in microfabrication techniques, penetration and ex- 
clusion of a single vortex by adjusting applied magnetic field 
has been observed repeatedly in many experiments lfTT3T4T2l . 
In Ref. [10 1, the penetration and exclusion of vortex is mea- 
sured by the Hall probe and in Refs. IfTTl Q~2) it is measured 
by a scanning tunnelling microscope (STM) tip. However, no 
statistics on the penetration and exclusion field has been car- 
ried out to investigate the thermal activation in these studies. 

For mesoscopic superconductors with several atomic lay- 
ers, the reduction of the dimensions promotes the thermal fluc- 
tuations. Even for mesoscopic superconductors made of con- 
ventional superconductors, G 2D = ^(k B TJ[H c (0)^(0) 2 d]) 2 * 
1CT 4 with the film thickness d ~ 1 nm, which is greatly 
enhanced by a factor (%/d) 2 » 1 compared to the bulk 
value fOl . By solving the Ginzburg-Landau equation and 
the Fokker-Planck equation analytically for a superconduct- 
ing disk, Pogosov found that the thermal activation of vor- 
tex over the BL barrier becomes feasible in mesoscopic 
superconductors! 13 1. The effects of thermal fluctuations on 
the penetration of single vortex have also been studied numer- 
ically by solving the time-dependent Ginzburg-Landau equa- 
tions with noise term| 14 1. They observed that the number of 
vortices in superconductors is fluctuating due to thermal acti- 
vation when G 3D is large in simulations. 

In this paper, we consider the thermally assisted penetration 
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FIG. 1. (color online), (a): A superconducting disk with radius of several coherence lengths and thickness of several nanometers. A STM 
tip is positioned at the center of the disk to check the presence of a vortex by measuring the zero-bias conductivity, (b): Profiles of the time 
dependent applied magnetic field. The maximal magnetic field is larger than the mean-field penetration field to initialize a vortex at the center 
of the disk. In the branch with increasing H, one measures the penetration field and in the branch with decreasing field, one measures the 
exclusion field, (c) and (d): Distribution of the (c) penetration field and (d) exclusion field P(H) at several temperatures. The mean-field 
penetration and exclusion field is represented by a dashed line. The distribution profiles are obtained using Eqs. (if), (28) and (30). 



and exclusion of a single vortex in a small superconducting 
disk. We find the distribution of penetration (exclusion) field 
when the applied magnetic field increases (decreases) adiabat- 
ically, as a function of disk size and temperature. The setup 
and main results are presented in Sec. II. In Sec. Ill, we de- 
rive the BL barrier for a superconducting disk by solving the 
London equation, while neglecting thermal fluctuations. From 
the BL barrier, we obtain the mean-field penetration and ex- 
clusion field for a single vortex. In Sec. IV, we calculate the 
thermal activation rate for vortex penetration and exclusion by 
solving the corresponding Fokker-Planck equation for vortex- 
as-a-particle and derive the corresponding distribution func- 
tions. The paper is closed by a discussion in Sec V. 



II. PROPOSED MEASUREMENT 

We consider a circular disk of several coherence lengths in 
lateral size and several nanometers in thickness d with d <K 
where the thermally assisted penetration and exclusion of a 
single vortex can be observed experimentally, see Fig. [JJa). 
An STM tip is positioned at the center of the disk to check 



whether the vortex is present or not, by measuring the zero- 
bias conductivity ifTTl . Without thermal fluctuations, the vor- 
tex penetrates into the disk when the applied field is larger 
than the mean-field penetration field. With thermal activation, 
the probability of vortex penetration is non-zero even when 
the applied magnetic is smaller than the mean-field penetra- 
tion field. To find the distribution of penetration field, period- 
ically varying magnetic fields are applied perpendicular to the 
disk. The setup with a static magnetic field has already been 
implemented in Refs. IfTTl fill [151 

For the vortex penetration, one changes the applied mag- 
netic field periodically and slowly, see Fig. [JJb). To present 
analytical calculations, we choose the sawtooth wave. One 
records the magnetic field when vortex enters in the branch 
with increasing field. By repeating the measurement, one then 
obtains a distribution of the penetration field at a given temper- 
ature. We calculate this distribution function analytically [see 



Eq. (30 1 below], which allows for a quantitative comparison. 
The measurement of the thermally assisted vortex exclusion 
is similar. It is determined from the branch with decreasing 
field, see Fig. [JJb). 
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III. ENERGY BARRIER 

In this section, we derive the mean-field surface barrier for 
vortex penetration and exclusion in a superconducting disk 
without accounting for thermal fluctuations in a static mag- 
netic field. The thickness d of the disk is much smaller than 
the London penetration depth d <k A. For such thin films, the 
screening of magnetic fields is weak and the effective penetra- 
tion depth is given by the Pearl length[16;l A — 2A 2 /d. The 
coherence length is also assumed to be much smaller than the 
penetration depth £ <sc A, and thus we can use the London ap- 
proximation, where the amplitude of the superconducting or- 
der parameter is uniform in space except for the vortex core. 
The magnetic field distribution h z (x,y) associated with a vor- 
tex is given by the London equation! 2 1 



K + — (V x g) • e 2 = <$ 6(r - r v ) , 

c 



(1) 



where g is the sheet current density, e z a unit vector along the 
Z direction and On = hc/2e is the quantum flux. Here r v is the 
position of the vortex and we only consider a single vortex. 
For the superconducting Pb film in Ref.[ 15 1 A as 20 yt/m. The 
lateral size of the film is much smaller than the Pearl length 
R <c A. In this case, the screening of the magnetic field is 
negligible. We can neglect the first term at the left-hand side 
of Eq. (JJJi. For convenience, we introduce a scalar stream 
function G, such that lfTTl 



Then we arrive at the total energy of a vortex in a circular disk 
of radius R 



U(r v ) 



8/r 2 A 



In 1 



n 

R 2 



+ ln 



R\ nR 2 H 



On 



1 

R 2 , 



(8) 

Equation ([8]) is valid when the vortex is not too close to the 
edge. For a vortex close to the edge, R - r v < one should 
account for the finite size of the vortex core. When R is of 
the order a different approach, such as direct solution of the 
Ginzburg -Landau equation, is needed lfl3l . The change of the 
magnetization when a vortex is at the center of the disk is 



M v = -d H U 



Opfl 2 
8ttA 



(9) 



Thus the penetration and exclusion of a single vortex can be 
determined by measuring M v using the Hall probe. ifTUl 

A vortex can be trapped in the disk if there is a local mini- 
mum at r v = 0. Local minimum develops when the magnetic 
field is larger than the mean-field exclusion field H e , 



Op 
nR 2 ' 



(10) 



Above H e there is a barrier for vortex exclusion. The position 
of the barrier is determined by the condition when dU /dp = 
with p = r v /R, which gives 



g = V x (Ge z ) . 



(2) 



At the boundary, the component of g normal to the bound- 
ary is zero, which requires that G is constant at the boundary. 
Without loss of generality, we take G — at the edge of the 
disk. Then Eq. ([TJ is reduced to the Poisson equation and the 
problem is equivalent to the one in electrostatics, 
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For a circular disk with radius R, we have lfT8l[T9l 



(3) 



2 , He 

It follows that the height of the potential at pi, is 
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(12) 



When U (ph) = 0, the barrier for vortex penetration vanishes, 
because the energy of a vortex at the edge is U (r v = R) = 0. 
This gives the penetration field H p 
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The energy of the vortex is then given by 
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(5) 



The energy e diverges at r = r v and to avoid the divergence, 
we introduce the standard cutoff £ near the vortex core, |r - 
r v | = Then the energy of vortex at r = r v reads 



e(r v ) 
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(6) 



We need also to add the interaction between the applied field 
H and the vortex, U z = -M V H with the magnetic moment of 
the vortex! 



M v 



dr(r x g) • e z . 
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R 

- —H e , 



(13) 



with e = 2.718. At fields above H p , nucleation of vortex into 
the disk becomes favorable. When the radius R increases, the 
surface barrier decreases thus H p also decreases. When U(p - 
0) = 0, a vortex at the center of the disk becomes the ground 
state, which gives another characteristic field H g 



nR 2 \4 



For R — 10£ = 300 nm, we have H e as 7 mT and H p 



(14) 



26 



mT. Several typical profiles of U(r v ) are displayed in Fig. [2] 
Slightly above H e , a local minimum appears at r v = 0, and 
U(R) < U(Q). When the field is increased above H g , the local 
minimum becomes the global minimum and U(R) > U(0). 
Above H p , the surface barrier for vortex penetration vanishes 
and the only stable solution is a vortex at the center of the disk. 
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FIG. 2. (color online). Profile of the energy barrier for a vortex in 
a circular disk with radius R/t; = 10 at several typical fields. The 
penetration field is H p = 0.074H c2 and the exclusion field is H e = 
0.02// t 2, with H C 2 = <t>o/(2n^ 2 ) the thermodynamic upper critical 
field. 



IV. THERMAL ACTIVATION 

When H approaches H p from below, the barrier decreases 
and thermal fluctuations will promote the penetration of a vor- 
tex. Thus in the presence of fluctuations, a vortex enters at 
magnetic fields below H p . For vortex exclusion, the vortex 
leaves the superconductor at fields above the mean-field value 
H e due to thermal activation. 

In the presence of time-dependent magnetic field as in Fig. 
[TJb), the energy barrier changes with time in this case. There- 
fore it is more appropriate to describe the vortex motion in 
terms of a Lagrangian. The Lagrangian of the vortex is 
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(15) 

The mass of the vortex is small and neglected. The dissipa- 
tion function is T) = i]r v 2 /2, where rj - <& 2 ) l(2nt; 2 c 2 'R n ) is the 
Bardeen-Stephen dissipation coefficient with % n the resistance 
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FIG. 3. (color online). Schematic view of thermally assisted vortex 
exclusion when H e < H < H g (a) and (b) penetration when H g < 
H < H p in a potential U(r v ). 



of the sample just above T c . r v < R is the coordinate of the 
vortex in the radial direction and r v is the radial vortex veloc- 
ity. The equation of motion for the vortex in the presence of 
thermal fluctuations is 



Wv 



d£(r v , t) 



+ F n , 



(16) 



where we have introduced a noise force F„ to account for the 
thermal fluctuations. F„ has a zero mean value (F„{t)) = 
and has a Gaussian correlator 

(F n (t, r v )F n ( f, r;,)) = 2rjk B T6(t - t')6(r v - r' v ), (17) 

where T is the temperature. When the change of the applied 
field is much slower than the thermalization time of the vortex 
(see Sec. V), we can use the adiabatic approximation to solve 
Eq. ( 16 1 and £,{r v , t) can be replaced by -U(r v ). A schematic 
view for thermally assisted vortex exclusion and penetration 
is depicted in Fig. [3] 

Here we have considered thermally assisted penetration and 
exclusion of a single vortex. For two vortices, the barrier is 
twice as large as that for a single vortex. Since the probability 
for thermal activation is much smaller than that of a single 
vortex, the simultaneous penetration and exclusion of multiple 
vortices can be avoided.! 



In Eq. ( 16 1, we have approximated the motion of the vortex 
in the disk by that of a particle in a one dimensional poten- 
tial. During the penetration or exclusion of vortex, the motion 
of single vortex is dictated by the force associated with the 
surface barrier. The diffusion of vortex along the azimuthal 
direction due to thermal fluctuations thus can be safely ne- 
glected. 



A. Thermally assisted exclusion 

Let us first calculate the vortex exclusion rate when H > H e 
using the Kramers equation for diffusion over barrier ll24ll 



i-U"(r m )U"(r h ) 
2jttj 



exp [-AU/k B T] 



(18) 



where {/"(rj) is the second derivative at the surface barrier 
and U"(r m ) is the second derivative at the energy minimum 
where the vortex resides initially. The barrier height is AU = 
U(rb) - U(r m ). The Kramers equation is valid when AU » 
k B T. 

Using Eq. 



18 i, we obtain the rate for the vortex exclusion 



7TT]R 2 



4 



2/(/-l) 2 exp 



(19) 



We 



with / = H/H e > 1, €a = 02/(8tt 2 A) and v = e c /(k B T). 
estimate e c ~ 2000 K for A w 20 jum. Hence v » 1 in the 
whole superconducting state T <T C « 7.2K. 

The characteristic frequency is to c = e c /{7U]R 2 ) as 10 GHz 
for R x 300 nm. The rate decreases when the lateral dimen- 
sion of the superconductor is reduced at a given H. More 
importantly, because e c and v decrease linearly with the film 
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thickness d, the thermal activation rate increases gigantically 
because v appears as an exponent. Note that the spurious de- 
crease of F e , when H is very close to H e , is an artifact, because 
Eq. ( 19 1 becomes invalid when AU ~ k B T . In reality, the rate 



r e increases monotonically when H approaches H e 



B. Thermally assisted penetration 

For vortex penetration, U"(r m ) is not defined at the edges. 
The Kramers equation becomes inapplicable, and an explicit 
solution of the Fokker-Planck equation is needed. 

We use the Fokker-Planck equation associated with Eq. 



dominant contribution is from the region near the barrier. We 
use the saddle point approximation 

^ * V{r b ) + \v"{r b ){r v - r b f, (25) 

with H"{rb) < 0. Furthermore we use the following approxi- 
mation because v » 1 



ex V [vW"(r b Xr v - r h ) 2 /2] 



Then the rate becomes 



In 



S(r v - r b ). (26) 



d r p(r v ,t) = -d r J(r v ,t), 



(20) 
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In 



(27) 



where p(r v , t) is the probability density of finding a vortex at 
r v and J is the probability current 

J(r v , t) = -Dexp[-U/(k B T)]d rr {exp[U/(k B T)]p}, (21) 

with the diffusion constant D = k B T '/rj. We use the standard 
boundary condition for the calculations of the rate [ 24 1, by as- 
suming a source of vortices at r,, = R and a sink at r v - 0. The 
total probability of finding a vortex in the region < r v < R 
is unity because vortex is initially thermalized at the edges, 



Let us compare Eq. (27i to the Kramers formula Eq. (18i. 
These two expressions coincide if U"(r m ) = 4nk B T '/£ 2 . For 
the vortex penetration, the vortex is thermalized at r v - R and 
gains kinetic energy k B T. Because of the uncertainty of the 
vortex position of the order one may assign an oscillation 
frequency of the order k B T/^ 2 . We then arrive at a consis- 
tent description based on the direct calculations of the Fokker- 
Planck equation and the Kramers equation. 

Finally, for the circular disk, we have the rate for vortex 
penetration 



pin 

> I d<f> I - 

Jo jR-$ 



r v dr v = 1, 
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p Rg\ R 



2V/C/-1) 



(28) 



which gives p(r v » R,t) = l/(2nR^). R-g appears in the lower with characteristic frequency io c - D/R£ * 10 GHz. Note 

that both V e and T p are strongly nonlinear functions of the 



bound of integral in Eq. ( 22 1 because of the uncertainty of the 
vortex coordinate due to the vortex core of finite size. We are 
interest in the probability for the vortex penetration, where 
initially no vortex is located at the center of the disk. To find 
the thermal activation rate, we can assume a sink at the center 
of the disk and use the boundary condition p(r v = 0, t) — 0. 
That means we restart the measurement once the vortex enters 
into the disk. 

Here we consider the regime where AU » k B T. In this 
case, the activation rate for vortex is small, thus J and d t p are 
small ll25l . Then the probability current is almost independent 
on r v , i.e., d r J = 0. Integrating Eq. (21 1 from r v = to r v — R 



and using the boundary condition, we have 



T 

Jo 



Dp(r v *R) = J exp[U/(k B T)]dr r , 



(23) 



where we have used U (r v = R) = because the vortex does 
not experience the surface barrier at the rim. Thus the thermal 
activation rate is given by 



radius R and the applied field H. 



C. Results 

To estimate the activation rate, we use the parameters from 
experiments 1151 . A — 150 nm, £ = 30 nm at 2 K and T c 
7.2 K. We take R = 300 nm, d = 2.5 nm and <R„ = 400 Q. 
The temperature dependences of ^(T) and A(T) are calculated 
with the BCS theory. The rate for the vortex penetration and 
exclusion is shown in Fig. |4] 

Experimentally, one measures the distribution of the pen- 
etration or exclusion fields when the magnetic field is swept 
periodically. The distribution of the penetration field can be 
obtained from the rate in the following way |f26|. The proba- 
bility for a vortex entering the disk after a waiting time t for 
the branch with increasing magnetic field in Fig. [TJb) is 



D 

F„ = 2nRJ = - 

p ^ 
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Jo 



W[H(t)] = 1 - exp 



exp[U/(k B T)]dr v 
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Jo 



r p [H(t')]dt' 



(29) 



(24) 



The characteristic energy scale of the barrier is e c . Since 
v = e c /(k B T) » 1, the integrand drops very rapidly and the 



Hence the probability distribution of the penetration field is 
dW T P [H] 



P(H) 



dH 



H 



■exp(-y), 



(30) 
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FIG. 4. (color online). Thermal activation rate for (a) the penetration 
of a vortex, (b) the exclusion of a vortex. The mean-field penetration 
and exclusion fields are also shown in the labels. Please note that 
the mean-field penetration field depends on temperature while the 
mean-field exclusion field is temperature independent. 



with 



y(H) 
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dH', 



(31) 



where H is the derivative of H with respect to time. For the 
penetration field measurement, H > 0. In a similar way, one 
obtains the distribution of the exclusion field for the branch 
with decreasing magnetic field in Fig. |TJb) 



P(H) = — : — exp 
\H\ 



Jh 



r e [H'] 

\H'\ 



dH' 



(32) 



For the exclusion field measurement, H < 0. Experimentally, 
the upper limit of the integral is replaced by the maximal mag- 
netic field used in experiment. The difference is negligible 
because T e ~ exp(-vH/H e ) w when H » H e . 

The distributions of the penetration field and exclusion field 
as a function of temperature are shown in Fig. [TJe) and (f). In 
the calculations, we use H - 14mT/s for penetration field 
and H — -14 mT/s for the exclusion field. As the temper- 
ature increases, the distribution profile becomes broad. The 
most probable penetration field decreases with temperature, 
and deviates from the mean-field value. As T increases, the 
exclusion field increases and also deviates strongly from the 
mean-field value. Therefore the hysteretic region between the 
vortex penetration and exclusion is reduced with increasing 
temperature. 



V. DISCUSSIONS 

Let us discuss the effect of time dependent applied field. 
First, the time dependent magnetic field induces an electric 
field and current, which generates heat. Secondly, the induced 



current exerts a Lorentz force on the vortex, thus it gives ad- 
ditional contribution to the surface potential. Thirdly, one 
should check the validity of the adiabatic approximation used 



in the evaluation of the rate in Eqs. ( 19 1 and (28 1 



The induced electric field according to the Faraday's Law 
is E(r) = -Hr/ilc), and the generated Joule heating is Q = 
cr q nH 2 R' [ dl(%c 1 ), where cr q is the conductivity of quasiparti- 
cles and is much smaller than the conductivity of the normal 
state. Using the heat balance equation, 



dC Pd d t T =%+a{T B -T) 



(33) 



the temperature increase after the system reaches a stationary 
state is AT = cr q H 2 R 2 d/{8c 2 a), where a is the heat exchange 
rate with substrate. Here C is the heat capacity and is the 
mass density of the superconductor. One can neglect heating 
if AT /T <K 1. For the sawtooth wave in Fig. [TJb) and (c), 
H is of order of a>fjH p with a> H being the period of the time- 
dependent field. To give an order of magnitude estimate, we 
take a « 10 7 W/m 2 ■ K and take the conductivity of Pb film in 
the normal state above T c , cr q d = 0.01 Q -1 . For u> H «: 1 THz, 
the heating effect is negligibly small. 

The induced current of Cooper pairs by the time-dependent 
magnetic field is accounted for by the last term on the right- 
hand side of Eq. (fSJ. The induced quasiparticles J q (r) = 
-cr q Hr/(2c) tilts additionally the potential for vortex pene- 
tration or exclusion. The change of the potential is 



AU„ = d I J a — dr = 
Jo c 



dcr q H<S>()R 2 
4c 2 



(34) 



It is required that AU q /e c <K 1. This condition is satisfied 
when o>h ^ 1 THz. 

To validate the adiabatic approximation employed in Eqs. 



( 19 1 and (28 1, the change of the applied field should be much 
slower than the thermalization time of vortex (or the relax- 
ation time of a vortex to reach the energy minimum), such 
that the vortex remains in equilibrium when the magnetic 
field changes ll27ll . From Eq. ( fl6] i, the thermalization time 
is J]/d 2 U ~ l/co c ~ 10~ 10 s for the parameters used in the pre- 
vious section. The adiabatic approximation thus is justified 
when a> H <s 10 GHz. 

Let us discuss the detailed form of the applied magnetic 
field. The rate description is valid when AU/(IcbT) » 1. 
The barrier near H p can be written as AU - e c (l - H/H p ). 
This gives an upper bound for the applied field 1 - H/H p » 
kBT/e c ~ 1/200 at T c . This means that the rate description is 
valid almost in the whole region of the magnetic field below 
H p . On the other hand, the period of the applied field should 
be large enough to count the event of vortex penetration in 
one period, which gives an upper limit for io H . The upper 
limit can be estimated as follows. The probability to observe 



vortex penetration according to Eq. (29 1 is, W = 1 - exp(-y). 



The rate can be written as F p « a> c exp(-v(l - H/H p )]. Then 



we have 



COrH,, 



y = — 



vH' 



exp 



-v 1- 



H 

H, 



VCU H 



exp(-10), (35) 
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where we have taken v(l - H/H p ) = 10, such that the rate de- 
scription is valid. In order to have large probability for vortex 
penetration in one period, y > 1, which gives lo h < 10 kHz. 
For such frequencies, the heating effect and tilt of the surface 
barrier are negligible. 

Some clues for the thermally assisted penetration and ex- 
clusion of a vortex in mesoscopic superconductors can be 
found by comparing two experiments IfTTI fl2l . In a recent ex- 
periment by Cren et. q/. lfT2ll . they measured the penetration 
and exclusion of the vortex in Pb films using an STM. No 
hysteresis in the penetration and exclusion of the vortex was 
observed. While in a similar experiment by Nishio et. al. IfTTI . 
hysteresis was clearly observed. These two experiments use 
Pb superconducting disks with similar thickness and radius. 
R a 2% « 60 nm in both experiments and d = 2.5 nm in 
Ref. IfTTI and d — 5.5 nm in Ref. [12]. The temperature in 
Ref. JT2) is 4.3K while in Ref.ETJ is 2K with T c * 7.2 K. 
According to Fig. [TJ(e) and (f), increasing the temperature 
reduces the hysteresis, which disappears eventually at higher 
temperature. The thermally assisted penetration and exclu- 
sion of vortex thus can qualitatively account for the difference 
in Ref. 1 121 and Ref. DTI. 

In summary, we have studied the thermally activated pen- 
etration and exclusion of a single vortex in a mesoscopic su- 



perconductor. We derived the energy barrier for a vortex as a 
particle using the London approximation, from which we ob- 
tained the mean-field penetration, exclusion field and thermo- 
dynamic critical field. We then calculated the thermal activa- 
tion rate for the vortex penetration and exclusion based on the 
Fokker-Planck equation and Kramers' escape rate in the adia- 
batic region of field change. Based on the activation rate, we 
obtained the distribution of the penetration and exclusion field 
due to thermal activation. Finally we proposed measurements 
of the distribution of penetration and exclusion magnetic fields 
to test our model of thermally assisted single-vortex motion in 
mesoscopic superconductors. 
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